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Q , Given a disk O in the plane called the objective, we want to find n small disks Pi, . . . , P n called 

■ the pupils such that U™j=i ^ ® Pj — 0> wnere © denotes the Minkowski difference operator, 

while minimizing the number of pupils, the sum of the radii or the total area of the pupils. This 
problem is motivated by the construction of very large telescopes from several smaller ones by 
so-called Optical Aperture Synthesis. In this paper, we provide exact, approximate and heuristic 
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Abstract 



a. 

■ solutions to several variations of the problem. 



1 Introduction 



The diameter of the pupil of a telescope is proportional to its resolution power. A simple calculus 
shows that we would need a telescope having a diameter of approximately 20m to observe the 
Earth from a high orbit [11]. Needless to say, such an instrument would not be adapted to the 
observation from space. In order not to build too large pupils, Optical Aperture Synthesis is 
adopted to synthesize (very) large pupils by interferometrically combining several smaller pupils [3] 
(see Fig. 1). The auto-correlation support (ACS) of a system of pupils denotes all the observable 
""^5 '. spatial frequency domain. 

The underlying problem can be stated in geometric terms as follows. Given an objective O 
supposed to be a disk, design a set of disks V = {Pi, . . . , P n } such that its ACS V covers entirely 
the objective while minimizing some cost function. Here T> = \Jiij = i(Pi Pj) where denotes the 
Minkowski difference operator. The cost function may include the number of pupils, the sum of the 
radii or the total area of the pupils, etc. This problem is a variant of the disk-covering problem. 
To the best of our knowledge, the variant we consider is new and the interferometry problem has 
not been considered before from a geometric perspective. This paper is a follow-up of our initial 
investigation [11]. The reader interested in the general disk-covering problem or some other variants 
can refer to [2, 6, 5]. 

The outline of this paper is as follows. In section 2, we introduce Apollonius diagrams (additively 
weighted Voronoi diagrams) which play a central role in our study, and use them to decide whether 
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Figure 1: Examples of using Optical Aperture Synthesis to synthesize large pupils [3] 




Figure 2: An Apollonius diagram of 8 disks in the Euclidean plane. The black disk has no cell. 



the objective is covered. Section 3 deals with the case of three pupils for which we provide an 
optimal solution. We describe in section 4 a constant-factor approximation algorithm for the case 
where the pupils are restricted to have the same radius. In section 5, we consider the centers of the 
pupils to be given and provide efficient algorithms to minimize the sum of the radii or the total area 
of the pupils under the constraint that the ACS covers the objective. Finally, section 6 considers 
the problem where the radii of the pupils are known but their positions are unknown. 

2 Apollonius diagrams and the decision problem 

2.1 Apollonius diagrams (aka Additively weighted Voronoi diagrams) 

Let V = {.Di, . . . , -Dat} be a set of N disks in the plane. We denote by q the center of Di and by pi 
its radius. Let ||.|| denote the Euclidean distance and dS denote the boundary of a subset of points 
S. The distance of a point x to the circle dD-i is defined as 

5i(x) = \\x - Ci\\ - pi. 

For a point x, 5i(x) is < 0, 0, > depending whether x lies inside, on the boundary of, or outside 
Di. The Apollonius cell of Di consists of the points whose distance to dDi is less than or equal to 
their distance to any other circle of T>: 

Ai = {x e R 2 | 5i(x) < 5j{x),j = 1, . . . , N}. 

Unlike the case of points, it is possible that a disk may have an empty cell. This happens 
when the disk is inside another disk. The one-dimensional connected sets of points that belong to 
exactly two Apollonius cells are called Apollonius edges, while points that belong to at least three 
Apollonius cells are called Apollonius vertices. The collection of the cells, edges and vertices forms 
the Apollonius diagram of V, denoted by Apo(V) (see Fig. 2). The Apollonius diagram Apo(V) 
can be computed in time 0(N log N) which is worst-case optimal [9], and robust and efficient 
implementations exist [1]. More information on Apollonius diagrams can be found in [4, 9]. We 
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Figure 3: A system of three pupils (left) and its ACS (right), the objective is represented by a thick circle. 



start by stating some properties of Apollonius diagrams. Let Bij define the bisector of two disks 
Di and Dj 

B i:j = {x eK 2 I Si(x) = Sj{x)}. 

Lemma 1. The restriction of <5j and 5j to B^ are unimodal functions. More precisely, these 
functions decrease linearly to a minimum and then increase linearly. 

Proof. Consider two disks -Dj and Dj with radii Pi,Pj and centers, w.l.o.g., Cj = (— c, 0) and Cj = 
(c, 0) . The bisector of Di and Dj is a sheet of the hyperbola whose equation is 

x^_ _ y 2 

9 9 9' 

where a = \pj — Pi\/2. Then the distance of a point with abscissa x on the hyperbola to Cj is a 
linear function of x: d = ±(ex + a), where e = - is the eccentricity of the hyperbola and sign ± is 
positive if pi < pj and negative otherwise. □ 

Corollary 2. Any arc pq contained in the edge of a cell Ai is included in the smallest disk of center 
Ci that contains p and q. 

Proof. Since the distance function to Di of the points on arc pq is unimodal by Lemma 1, it reaches 
a maximum at p or q. Hence any disk with center Cj that contains p and q covers the whole arc. □ 

Corollary 3. The Apollonius cell Ai is included in the disk centered at Ci that contains the set of 
its vertices. 

Proof. If Ai is unbounded we are done. Otherwise, as Ai is star-shaped [4], it is included in a disk 
if its edges are. Applying Corollary 2 to all edges of Ai concludes our proof. □ 

Let 5x>(x) denote the smallest distance of x to the disks of V, i.e., 5x>(x) < 5i(x) for any 1 < i < N 
and equality holds iff x E Ai. We see that Sx>(x) < when x lies inside the union of the disks of T>. 

2.2 The decision problem 

Let V = {Pi, . . . , P n } be a set of n disks called the pupils and O be a disk of radius R centered at 
the origin called the objective. The ACS of V is T> = \J i j=\(Pi &Pj)- The decision problem consists 
in determining whether O is covered by T>. 

Let Ci and pi denote the center and the radius of pupil Pi and let Dij = Pj Pj . It is not difficult 
to see that D^ is a disk with center Cy = Q — Cj and radius pij = Pi + Pj. Moreover, V = U" J=1 Dij 
(see Fig. 3). 

If the radius pi of some pupil Pi is greater than half the objective's radius R, Da covers O. We 
assume in the sequel that the pupils all have a radius at most ^ which implies that all disks of T> 
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Figure 4: Left: A set of three pupils whose ACS does not cover the objective. The x-marks correspond the 
vertices of Vij of which some lie outside the union of disks. Right: The set of pupils with the same position 
but radii enlarged by a* as computed by Algorithm 1. All vertices of Vij are inside V and the objective is 
covered. 



have radii smaller than R. We write Ay for the cell of JDy in the Apollonius diagram of T>. Let Vy 
denote the set of vertices of Ay inside O and the intersection points of dAij with dO. We denote 
by N = n 2 the number of the disks of T>. It can be argued that the cardinality of all Vy is O(N). 
The following shows a necessary and sufficient condition for covering O by Ui?=i^y ( see Fig. 4). 

Lemma 4. O C T> iff Vy C Dy for all i, j = 1, . . . , n. 

Proof. First we argue that O C T> iff Ay HOC £)y for all i,j = 1, . . . ,n. Since the set of Ay- 
forms a decomposition of the plane, Ay HOC £)y, i, j = 1, . . . , n, implies that O C U" J=1 Dy = 2?. 
Conversely, suppose that O C D and p € Ay n O, we will show that p E .Dy. Indeed, peOCP 
implies < 0. Together with p G Ay, we conclude <%(p) = < which implies that p is 

inside -Dy. 

We show next that Ay ROC Dy is equivalent to Vy C Z?y by proving that a disk A centered 
at cy covering Vy covers also Ay n O. We first observe that the edges of Ay with both endpoints 
in O are covered by A by Corollary 2. It remains to verify that the intersection points of <9Ay 
with dO and the arcs linking them are also in O. Consider two such points p and q consecutive 
along the boundary of O. Call pip2 and q±q2 the two Apollonius edges that intersect dO at p and q 
respectively. Suppose pi,qi G O and P2,q2 $ O, which implies that pi,p,q,qi belong to Vy. Since 
pi and p lie on edge P1P2, and q and q± are contained in q±q2, A will cover the arcs p\p and qq\ by 
Corollary 2. It thus remains to show that the circular arc pq of O is included in A, which is true 
since p,q G Dy whose radius has been assumed to be smaller than the radius of O. □ 

The following simple result is important in sections 5 and 6. 

Corollary 5. Given a configuration of pupils with the corresponding sets Dij and Vij. We move/resize 
the pupils such that each new disk includes Vij . Then, O is covered by Uij=i -^y ■ 

Proof. Since Vij C is equivalent to Ay HOC (see the proof of Lemma 4) and the sets 
Ay n O cover O, the objective is covered by [J™j =1 D'-. □ 

Lemma 4 gives us a simple 0(AHog iV)-time algorithm that solves the decision problem. It still 
works when we replace Apollonius diagrams by power diagrams. The reason of using the formers 
will be seen in section 5. 

3 Problem with three pupils 

A configuration of pupils is called valid if its ACS covers the objective. In this section, we want 
to minimize the sum pi + P2 + P3 among the valid configurations. Let denote by I the line passing 
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Figure 5: A configuration of three pupils and the upper part of its ACS 



through C23 and C32. Since the disks and the objective are symmetric about the origin, it suffices to 
consider only one half-plane bounded by /. 

Lemma 6. Among the valid configurations, those in which one radius is half of the objective's 
radius R and the other two are zero are optimal. 

Proof. It is straightforward to see that such configurations are valid. Consider now a configuration 
in which p\ + p2 + p3 < R/2. We will prove that it cannot be a valid configuration. Indeed, suppose 
w.l.o.g. P\ has the largest radius among three pupils. Then D\\ is the largest disk among Dn,D22 
and -D33 and its radius 2p\ is smaller than R. Let pi,qi,q2,P2 be the intersection points from left 
to right of dO and dDn with I (see Fig. 5). If segment pTqi is covered by D23, then the diameter 
of -D23 is at least the length of pTqi, i.e., 2(p2 + ps) > R — 2p± which implies p\ + p2 + pz > R/2 
(a contradiction). The case where p~2q~2 is covered by D23 is symmetrical. We can therefore assume 
that -D23 does not cover p±qi nor p~2q~2, and, by symmetry, the same holds for D^2- Without loss of 
generality, we can assume that D12 contains p\ or q\ and that D13 contains P2 or 52 • We denote by c 
the midpoint of the arc p\p2 of dO. The distance of c to pi, qi,P2, <?2 is at least \J R 2 + (2pi) 2 > R. 
Then c is not included in neither D\2 nor D13 whose diameters are smaller than R. It is not included 
in D23 and .D32 either since the distance from c to C23 and C32 is at least R and the radii of D23 and 
D32 are less than R. Hence, the configuration is not valid. □ 

It is interesting to see from the above lemma that configurations of three pupils consisting of a 
pupil of radius R/2 and two points are optimal, whatever the position of the pupils may be. 

4 An 8^— approximation to the smallest number of the pupils of 
the same radius 

In this section, we restrict to the case p\ = . . . = p n = p/2, then the disks Dij have the same radius 
p. We want to find an upper bound for n to cover an objective of radius R. As the number of disks 
is n 2 , a lower bound \R/p] is easily obtained. 

Let p be any prime number, we start by stating a basic property of p 

Fact 7. Let k,l G Z such that gcd(p,k) = 1, there exists an integer < i < p satisfying ik = I 
(mod p). 

Theorem 8. {xj — Xj \ i, j = 0, . . . ,4p — 1} D {x G Z, \x\ < p 2 } where 



kp + (- 

Xk +P, 



fc(fc + l) 
2 



mod p) 



Xk+2p 



for k = 0...,2p-l. 
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Proof. Let x be an arbitrary integer between and p 2 — 1, then x can be written as kp + l for some 
< k, I < p. Let Xi = Xk+i — Xi for i = 0, . . . ,p — 1, we observe that 

(* - l)p < Xi < (k + l)p. (1) 

Xj = Xo + ik (mod p) 

By Fact 7 there exists some < i < p such that Xi = I (mod p). Hence together with (1) the 
difference of either x^+i or x^+i+2p with will be x. The only case where Fact 7 does not apply is 
when k = 0. In this case choose fc = 1 instead and easily see that the set {xj+i — X{ + 2 P } U {xj+i+2p — 
Xj+2p} generates all integers 1, ... ,p — 1 and hence contains x. □ 

The above set should not be confused with Golomb ruler [7] and the set defined by Erdos and 
Turan [8] since in the latter sets, the differences between any pair of distinct elements must be 
unique but do not generally cover all points 1, . . . ,p . 

Suppose, w.l.o.g., radius of the disks p = ^= and R = p 2 for some prime p. Let S = {x G 

^ 2 I | \x\ |oo < P 2 }- We see that the disks of radius whose centers cover S are sufficient to cover 
completely the objective. In other words, we want to find n centers of pupils q € I? such that 

{cj — Cj | 1 < i, j < n} 5 <S 

Corollary 9. \8y/2R/ p] pupils of radius p are sufficient to cover an objective of radius R. 

Proof. The set of pupils is constructed as follows: q = (xij_\,Xi mo d 4p) for i = 0, . . . , lQp 2 — 1. By 

applying Theorem 8 first for x-coordinate and then for y-coordinate, we see that these 16p 2 pupils 
are able to cover any element of S thus the objective of radius R. As R = p 2 and p = we yield 
the upper bound. □ 

The following is an immediate consequence of Corollary 9 and the lower bound observed earlier. 

Corollary 10. There is an 8y/2— approximation algorithm to cover the objective of radius with the 
smallest number of pupils of the same radius. 



5 The fixed-center problem 

In sections 5.1 and 5.2, the centers of the pupils are fixed and we present two heuristic algorithms 
for optimizing the radii among the valid configurations. Both algorithms are based on the fact that 
the circle of center Cy and radius pij + 5ij(p) passes through p. Then we provide an approximation 
algorithm with a given error bound and compare it with the heuristic algorithms. We end up the 
section with a method to maximize the objective while keeping fixed the radii as well as the positions 
of the pupils. 



5.1 A simple optimization problem 

If we increase each of the radii of the pupils by a real number a/2, the radii of the disks Dij then 
increase by a and Apo(T>) remains unchanged. Hence there exists a minimal value of a, denoted 
a*, for which the objective is covered by the union of the new (enlarged) disks. 

The following shows that a* can be computed exactly in 0(N log N) time (see Algorithm 1). 
We recall that Vij is the set of vertices of Aij inside O and the intersection points of dAy with 80. 

Lemma 11. 

a* = max max Sij (p) 

Proof. It is easy to see that maxjj max pe y i;) 5ij(p) is the minimal value of a for which Vij C Dij. 
The result follows from Lemma 4. □ 
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Algorithm 1 Compute a* 

1: a* < oo 

2: compute Apo(V) and Vij 

3: for all cells A^ of Apo(V) do 

4: for all x G V^- do 

5: a* <— max(a*, <%(o;)) 

6: end for 

7: end for 

8: return a* 



5.2 Minimizing the sum of the radii of the pupils 

We consider now the more difficult problem of optimizing the sum of the radii of the pupils and 
propose a heuristic solution that turns out to perform well in practice. 

Instead of increasing the radii of the Pj by a same amount as in the previous subsection, we 
consider them as n variables. Algorithm 2 below proceeds in two main steps. First, we compute 
minimal quantities, denoted a^j, by which the radii of the Dij must be enlarged/reduced so as to 
satisfy Lemma 4 (lines 3-9). This step is similar to Algorithm 1. Thanks to the fact that the 
already visited a^ necessarily increase, the initial Vij will be covered upon termination by the disks 
(which are augmented by aij). The objective is then covered by \Ji j=i^ij according to 
Corollary 5. Finally, we want to minimize the sum of the radii of the P* under the constraint that 
p* + p* must be at least the radius of D- • (line 10): 



min p i 

i=i 

s.t. p* + p* > {pi + pj) +a>ij, i,j = l,...,n (*) 

p* > 0, i = l,...,n. 

Here, pi are the radii of the initial pupils Pi and hence known. This is a linear program whose 
feasible set is non-empty and bounded. Thus, there exists an optimal solution. 



Algorithm 2 Minimize the sum of the radii of the pupils 



e <— any small positive constant 
repeat 

< oo, i, j = 1, . . . ,n 

compute Apo(V) and Vij 
for all cells A^ of Apo(V) do 
for all x G Vij do 

a^ <- max(a ij ,S ij (x)) 
end for 
end for 

compute {p*}j=i,..., n by solving the linear program (*) 

i=l Pi ~ Ei=l Pi 
Pi^ p*, i = !,■■■, n 

until err < e except for the first iteration 

return {p*}j=i,..., n 



Note that we need to update the Apollonius diagram since the pupils' radii change after each 
iteration of the repeat loop. 

Lemma 12. Algorithm 2 always terminates. 
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Proof. The initial Vij is included in by the construction of a^. According to Corollary 5, O 
is therefore covered by [j - Dij after the first iteration. Hence, we may assume that the objective 
is covered. In this case, Lemma 11 implies that no a^- is positive which shows that, at each step, 
P* + Pj < Pi + Pj and hence Yli=i Pi — Y17=i Pi- Since Yli=i Pi ^ s positive, Algorithm 2 necessarily 
terminates after a finite number of iterations. □ 



Minimizing the total area of the pupils: Replacing the objective function Yli=i Pi m (*) with 
7T J27=i Pi 2 yields a quadratic program which minimizes the total area of the pupils. 

Additional constraints: In addition to covering the objective, we can also bound the radii of 
the pupils and forbid any overlap among the pupils. This can be done by adding the following 
constraints to the linear program (*) 

P* +P* < \\ci-Cj\\, l<i<j<n, 
minjradius < pi < max -radius, i = 1, . . . , n. 

Algorithm 2 has been implemented and appears to work well in practice. Fig. 6 compares the 
results of Algorithms 1, 2 with the optimal solution computed by the following exhaustive search 
method. 

Exhaustive search algorithm: If the radii of the pupils are assumed to be integer multiples of 
a small number 9, then the exhaustive search methods can be applied and the optimal solution in 
the continuous case must be at least the solution found by these methods minus n9. We hence have 
an approximation algorithm within a given error bound. 

5.3 Maximizing the objective 

Now we keep the pupils fixed (radii and positions) and maximize the radius of the objective under 
the constraint that it is covered by the union of the disks. 

Proposition 13. // an edge pq of Aij cuts dDij at a point x / p and q, then there is a point x' on 
pq that is close to x and not contained in T>. 

Proof. From the fact that %(.) is a unimodal function and 5ij{x) = Sx>(x) = 0. □ 

The following corollary, whose proof is referred to the full version of the paper, computes the 
maximal radius R* of the objective for which it is covered by V. 

Corollary 14. If Da C An for some i = 1, . . . , n then R* = 2p, L . Otherwise, 

R* = min min ||x|| 

ij xGdAijddDij 

6 The fixed-radius problem 

In this section we fix the radii and propose a heuristic algorithm for moving the set of pupils so that 
its ACS covers the objective. Our algorithm is based on Corollary 5. More precisely, we want to 
capture the point sets Vij by the disks D^. Given a set of points P and a disk, the optimal center 
position for the disk to cover P is the point that minimizes the maximal distance to any point of P 

min max \\x — p\\. (2) 

peP 

This is the so-called smallest enclosing disk problem and a linear algorithm to compute exactly the 
disk center can be found in [10]. Unfortunately, function (2) being non-differentiable makes it hard 
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Figure 6: Initial configuration of 5 pupils (a). Results after applying Algorithm 1 (b), Algorithm 2 (c) and 
the exhaustive search algorithm (d). The total areas of the pupils in (a), (b), (c) are 35.6571, 27.1062 and 
19.8421 respectively. The optimal solution must be at least 19.572 as computed by the exhaustive algorithm. 
The area of pupils in (e) is only 16.4793 when we move the pupils by the algorithm in section 6 and then 
apply Algorithm 2. 




Figure 7: a) The pupil centers are initially placed about a planar line, b) Dotted curves illustrate the 
movements of the pupils after iterating 24 times the algorithm in section 6 when the union of disks cover 
completely the objective, c) A configuration of 9 pupils d) Result obtained by iterating 9 times the algorithm. 




Figure 8: a) A configuration of 10 pupils, b) Result after applying Algorithm 2. c) Result after moving the 
pupils by the algorithm in section 6 and then applying Algorithm 2. The total areas of the pupils in (a), (b) 
and (c) are 57.0827, 48.0171 and 26.8088 respectively. 
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to apply to our problem. Another approach is to minimize the sum of the squared distance from 
the center to each point of P 

min — p\\ 2 . 

This function is convex and attains its minimum at the barycenter of P. Our algorithm works 
as follows. We begin with a given configuration of pupils, compute the set Vij and move the pupils 
to minimize the following function 

n 

mmE £ \\(c* - c*) - pf 

Here the centers c* of the pupils are variables and we recall that c* — c* becomes the center of disk 
D*y The objective function being the sum of convex functions, is thus convex. We can update the 
sets Vij and iterate the algorithm until we obtain the desired result. As shown in Fig. 7, the initial 
configuration is not critical. The algorithm can also be used as a preprocessing step to improve 
Algorithm 2 (see Figs. 6e and 8). 

ACKNOWLEDGMENT. We thank Helmut Alt, Giinter Rote and Mariette Yvinec for helpful 
discussions and careful proofreading of early drafts of this paper. 



References 

[1] CGAL release 3.2.1, http://www.cgal.org/. 

[2] H. Alt, E.M. Arkin, H. Bronnimann, J. Erickson, S.P. Fekete, C. Knauer, J. Lenchner, J.S.B. 
Mitchell, and K. Whittlesey. Minimum-cost coverage of point sets by disks. In Proc. Symposium 
on Computational Geometry, pages 449-458, 2006. 

[3] P. Blanc, F. Falzon, and E. Thomas. A new concept of synthetic aperture instrument for high 
resolution earth observation from high orbits. In Disruption in Space, 2005. 

[4] J-D. Boissonnat, C. Wormser, and M. Yvinec. Curved Voronoi diagrams. In J-D. Boissonnat 
and M. Teillaud Eds., editors, Effective Computational Geometry for Curves and Surfaces, 
chapter 1. Springer, 2006. 

[5] L. Booth, J. Bruck, M. Franceschetti, and R. Meester. Covering algorithms, continuum per- 
colation and the geometry of wireless networks. Annals of Applied Probability, 13(2):722-741, 
2003. 

[6] J. Cortes and F. Bullo. Coordination and geometric optimization via distributed dynamical 
systems. SI AM Journal on Control and Optimization, 44:1543-1574, 2005. 

[7] A. Dimitromanolakis. Analysis of the Golomb ruler and the Sidon set problems, and determi- 
nation of large, near-optimal Golomb rulers. Master's thesis, 2002. 

[8] P. Erdos and P. Turan. On a problem of Sidon in additive number theory, and on some related 
problems. Proc. London Math. Soc, 16:212-215, 1941. 

[9] M. Karavelas and M. Yvinec. Dynamic additively weighted Voronoi diagrams in 2d. In Proc. 
10th European Symposium on Algorithms, pages 586-598, 2002. 

[10] M. Overmars M. de Berg, M. van Kreveld and O. Schwarzkopf. Computational Geometry: 
Algorithms and Applications. Springer- Verlag, 2nd edition, 2000. 

[11] T. Nguyen, J-D. Boissonnat, P. Blanc, F. Falzon, and E. Thomas. Pupil configuration for 
extended source imaging with optical interferometry: A computational geometry approach. In 
Proc. IEEE Int. Conf. on Acoust., Speech and Sig. Proc, volume 2, pages 793-796, 2006. 



10 



